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Large-area single-layer WS2 is grown epitaxially on Au(lll) using evaporation of W atoms in a 
low pressure H2S atmosphere. It is characterized by means of scanning tunneling microscopy, low- 
energy electron diffraction and core-level spectroscopy. Its electronic band structure is determined 
by angle-resolved photoemission spectroscopy. The valence band maximum at K is found to be 
significantly higher than at T. The observed dispersion around K is in good agreement with density 
functional theory calculations for a free-standing monolayer, whereas the bands at T are found to 
be hybridized with states originating from the An substrate. Strong spin-orbit coupling leads to a 
large spin-splitting of the bands in the neighborhood of the K points, with a maximum splitting of 
419(11) meV. The valence band dispersion around K is found to be highly anisotropic with spin- 
branch dependent effective hole masses of 0.40(02)me and 0.57(09)me for the upper and lower split 
valence band, respectively. The large size of the spin-splitting and the low effective mass of the 
valence band maximum make single-layer WS2 a promising alternative to the widely studied M0S2 
for applications in electronics, spintronics and valleytronics. 

PACS numbers: 73.22.-f, 73.20.At, 79.60.-i 


I. INTRODUCTION 

The discovery of graphen^i® established the possibil¬ 
ity to obtain stable two-dimensional solids, and it was 
soon realized that layered materials other than graphite 
can be used as bulk parents for novel two-dimensional 
materialJ^l. As in graphite, the weak van der Waals in¬ 
teractions between adjacent atomic layers in alternative 
parent materials, such as transition metal dichalogenides 
(TMDCs), permit mechanical exfoliation of a single-layer 
(SL) with electronic properties that differ in subtle but 
important ways from those of the bulk material. One 
example is the indirect to direct band gap transition ob¬ 
served in case of MoS and other TMDC^. 

SL TMDCs consist of a layer of metal sandwiched 
between two layers of chalcogens. These materials ex¬ 
hibit a variety of electronic properties, ranging from in¬ 
sulating to metallicP. Probably the most studied SL 
TMDCs are molybdenum dichalcogenides, which are 
attractive mat erials for electronic applications such as 
transistor j^ l ^^ l, diodephotoemitting device^^, solar 
cell^I^ and memristor^I^. Furthermore, the unique spin 
texture of both conduction and valence band makes 
SL TMDCs well-suited for studying quantum degrees 
of freedom such as spin or valley pseudospin or their 
interaction^^^JhZI, In the case of the tungsten dichalco¬ 
genides, much stronger spin-orbit coupling is expected 
than in the case of the Mo-based analogues, and the 
properties just enumerated should thus be more stable 
at room temperaturd^^^^. Naively, a strong spin-orbit 
splitting of the bands can be expected to result in an in¬ 
creased band curvature near the top of the valence band 
and hence to a reduced effective mass; and indeed, WS 2 
is predicted to be the best material among all of the 
TMDCs for a transistor channel, due to its low effective 
hole mas^^^^f^. 


Many proposed approaches to characterizing the elec¬ 
tronic properties of SL TMDCs, as well as many potential 
applications, require large area and high quality samples. 
In the present work, we introduce a method for the epi¬ 
taxial growth of WS 2 on Au(lll). We study the growth 
and structure by scanning tunneling microscopy (STM), 
low energy electron diffraction (LEED) and core level 
spectroscopy. The high quality of the SL WS 2 obtained in 
this procedure permits an investigation of the electronic 
structure by angle-resolved photoemission spectroscopy 
(ARPES). The results confirm the expected strong split¬ 
ting of the SL WS 2 valence band at K, indications of 
which had already been seen in early experiments on sin¬ 
gle layer^^. The splitting is found to be almost three 
times larger than the spin-splitting determined for SL 
MoSj^with a strong warping of the constant energy con¬ 
tours around the K point. The expected lower effective 
mass near the valence band maximum is also confirmed. 


II. METHODS 

Growth and measurements were performed at the 
SGM3 endstation at the ASTRID2 synchrotron radiation 
facilit}ESl. Prior to growth, the Au(lll) single crystal 
substrate was thoroughly cleaned by repeated cycles of 
Ne ion sputtering (E=0.75 keV) and annealing (600°C) in 
ultra high vacuum (UHV) until the regular herringbone 
reconstruction was observed by STM. In the first step 
of the growth procedure, tungsten was evaporated onto 
the clean surface using a commercial e-beam evaporator 
charged with a 99.95% purity W rod. During the evapo¬ 
ration, the sample was exposed to 99.6% purity H 2 S us¬ 
ing a homemade nozzle situated ~ 1 cm from the sample. 
This allowed for a local pressure at the sample face that 
was higher than the background pressure in the chamber 
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(which was maintained at ~ 10“^ mbar). In the second 
step, the sample was annealed to 650°C for 30 minutes, 
maintaining the H 2 S atmosphere. This procedure is sim¬ 
ilar to that used in Ref. [27|to obtain WS 2 nano clusters. 
In order to obtain a large coverage of WS 2 , however, it 
is necessary to repeat the above sequence. The num¬ 
ber of growth cycles determines the coverage of the SL 
WS 2 film. A last anneal at 750° C under UHV conditions 
was used to desorb any residual contamination from the 
growth process. The sample quality was examined by 
STM, LEED and core level spectroscopy using a photon 
energy of 140 eV. The reported core level binding energies 
have been calibrated using the bulk component binding 
energy of the clean Au(lll) substrat^^. The character¬ 
ization of the film was performed at room temperature. 

ARPES measurements were carried out at a temper¬ 
ature of ~ 110 K, with energy and angular resolutions 
better than 20 meV and 0.2°, respectively. Measurements 
were taken with photon energies in the range of 14-80 eV. 
The WS 2 -related bands did not show a kz dispersion in 
this energy range, confirming the two-dimensional char¬ 
acter of the system (photon energy scans not presented 
here). 


III. RESULTS AND DISCUSSION 

Representative STM images of SL WS 2 on Au(lll) 
are shown in Eig. The small lattice mismatch be¬ 
tween the substrate (a^^i(iii) = 2.88 A) and SL WS 2 
{ciwS 2 — 3.15 A)P3 leads to a modulation of the local 
densi^ of states, resulting in the observed moire pattern 
(Eig. [^b)). The apparent height of the WS 2 (~ 3.1 A) 
and the moire superstructure periodicity (^ 31 A) are 
determined from the respective line profiles in Eig. [^c) 
and (d), and are c onsist ent with the results for WS 2 and 
M 0 S 2 nano-islandPI^. The SL WS 2 coverage of the 
sample shown in Eig. [^is approximately 0.7 monolayers 
(ML). At this coverage, some bilayer regions begin to be 
observed. This amount of bilayer coverage is not large 
enough to be detectable in ARPES, where the formation 
of a bonding / anti-bonding splitting of the topmost va¬ 
lence band leads to a distinct difference between the elec¬ 
tronic structure of bilayer and SL TMDC^ (see further 
discussions of this below). 

The formation of the moire pattern can be followed by 
LEED measurements taken at lower coverages as shown 
in Eig. Eor a coverage of ~ 0.15 ML, two concen¬ 
tric hexagonal patterns dominate (Eig. [^a)). Based 
on length of the unit cell vectors, the outer and in¬ 
ner hexagon can be assigned to Au(lll) and SL WS 2 , 
respectively. The measured reciprocal unit cell vector 
ratio excellent agreement with 

the crystal unit cell vector ratio 

the case of the higher coverage (~ 0.3 ML) shown in 
Fig-ib), the moire pattern is clearly observed in addi¬ 
tion to the spots caused by the separate reciprocal lattice 
vectors of Au(lll) and WS 2 . The presence of the moire 
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FIG. 1: (color online) (a) Large-area STM image, showing 
regions of bare Au(lll), SL WS 2 (coverage of ~0.7 ML) and 
bilayer WS 2 . (b) Small-area image emphasizing the moire 
structure which has a measured periodicity of ~ 31 A. (c) 
Apparent height profile along line 1 in panel (a), showing an 
approximate height of 3.1 A for SL WS 2 . (d) Line profile 
acquired along line 2 in panel (b) showing the moire period¬ 
icity. Scanning parameters : (0.06 V, 0.2 nA) and (0.62 V, 
0.28 nA), in (a) and (b) respectively. The STM images were 
analyzed using the free WSxM softward^. 


pattern is consistent with the STM results reported in 
Eig. [^b). Its absence for the lowest coverage can be 
explained by the fact that the growth starts with the for¬ 
mation of small islandP^, which are initially too small 
to establish a moire pattern. Similar LEED patterns 
have also been observed for graphene grown on transi¬ 
tion metal^^. 

Fig-i shows angle-integrated core level spectra taken 
in order to follow the local chemical changes during the 
growth of SL WS 2 /Au( 111 ). All spectra were analyzed 
by the subtraction of a Tougaard-type ba ckg rouncP^ and 
subsequently fitted by Lorentzian profiled. Before the 
growth, the An 4f core levels of the clean Au(lll) crystal 
were measured as a reference (Eig. [^a)). A fit reveals a 
bulk component and a surface-shifted component at lower 
binding energy, in good agreement with the literatur^^. 
The growth of approximately 0.7 ML of SL WS 2 was 
found to affect the Au 4f core levels as shown in Eig. 
ib). The surface component intensity is largely reduced 
and shifted towards the bulk component by 58(32) meV 
with respect to clean Au(lll). It is difficult to analyze 
this spectral change in detail, as there are several factors 
contributing to it. The first is the actual suppression or 
shift of the surface core level component due to the in¬ 
teraction with the sulphur. The second is the remaining 
presence of clean Au(lll) patches that should give rise to 
a much weaker but largely un-shifted surface component. 
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FIG. 2: (color online) FEED images at kinetic energy of 
118 eV for samples of different coverages, (a) Eor a low cov¬ 
erage sample 0.15 ML) two concentric hexagonal patterns 
are observed originating from the substrate and SL WS 2 . The 
reciprocal unit cell vectors are indicated with arrows, (b) In 
case of samples with higher coverage (~ 0.3 ML) a moire pat¬ 
tern is observed due to sufficiently large island size. The inset 
presents a magnified moire pattern observed around one of 
the main spots. 


The W 4f core level spectrum consists of a spin-orbit 
split doublet with each peak showing an intense high 
binding energy component and a weak low energy compo¬ 
nent. The observed binding energies for the main compo¬ 
nents are E 5/2 = 34.85(02) eV and E 7/2 = 32.72(01) eV. 
The spin-orbit splitting of the states is thus 2.13(02) eV. 
These values are con sistent with measurements on other 
WS 2 systernP^mm^ However, there is a clear difference 
between the spectra here and those previously reported 
for smaller nano-scale WS 2 islands on Au(lll)P^, for 
which two almost equally strong components were found 
in each spin-orbit split peak of the W 4f core level. Those 
findings were interpreted in terms of core level compo¬ 
nents from the edges as well as from the basal plane of 
the WS 2 nano clusters, consistent with findings for M 0 S 2 
cluster^^. This interpretation is supported by our result 
of only one main component, as the number of edge atoms 
with lower coordination is small for the high coverage re¬ 
alized here. The weak low binding energy component 
observed in Fig. |^c) might stem from edge atoms, or 
from the presence of W atoms in lower oxidation s tates, 
probably due to not fully sulfidized regions 
The presence of metallic W(0) or W(+II) is not observed. 

Fig. 0a) shows the electronic structure of SL 
WS 2 /Au( 111 ) along different high symmetry directions 
of the Brillouin zone as observed by ARPES. Correspond¬ 
ing constant energy contours are presented in Fig. [^b). 
Sharp features attributed to Au(lll), such as the surface 
stated, sp bulk bands and projected bulk gap edgeP^, 
are marked in Fig. [^a). 

The uppermost valence band of SL WS 2 is discernible 
for binding energies larger than 1.2 eV, especially be¬ 
tween F and K and towards M. The band is signifi¬ 
cantly sharper within the projected bulk band gaps of 
Au around K (outlined by dashed lines) than near F 
and M. Near M the WS 2 states strongly hybridize with 
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FIG. 3: (color online) (a-b) Au 4 / core levels for clean Au(lll) 
and J^O.TML of SL WS 2 /Au(lll), respectively. Blue and 
red peaks correspond to bulk and surface components, re¬ 
spectively. (c) W4/ core levels from SL WS 2 . Small peaks 
observed at lower binding energy (in magenta) indicate the 
presence of W atoms with lower oxidation states (WS 2 -cc)- 


the gold continuum and cannot be discerned. The max¬ 
imum of the valence band is found to be situated at K 
(^ 0.51 eV higher than at F), consistent with the ex¬ 
pected direct band gap at K. The observation of a single 
valence band at F with a higher binding energy than at 
K also rules out a significant contribution from bilayer 
WS 2 : in bilayer TMDC systems, the valence band near F 
shows a bonding / anti-bonding splitting that would be 
observable W ARPES as two bands, in contrast to what 
is seen her^. 

The overall observed dispersion of the band structure 
is in good agreement with calculations for a free-standing 
SL WS 2 (red dashed lines)P^, notably around K. Diver¬ 
gence from these calculations in terms of a shift towards 
higher binding energy is seen near F, similar to what has 
previously been observed for SL MoS 2 /Au(lll), where 
it was attributed to an interaction with the substrat^^. 
The bands around F are a mixture of out-of-plane W 
d ^2 and S Pz orbitals, which can participate in bond for¬ 
mations with the Au and as a result change the band 
dispersiorP^. In contrast to this, the valence bands at 
the K points are derived mainly from in-plane W 
and <1^2_y 2 orbitals. These are not only less likely to 
be affected by the interaction with the substrate, but 
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the states nearby are also situated in a projected bulk 
band gap. Note that this situation is distinctly different 
from the case of graphene where the bands near K are 
formed from out-of plane tt orbitals. This difference in 
orbital character near the K point also explains the ab¬ 
sence of replicas induced by the moire superstructure for 
SL WS 2 /Au( 111 ), something that is observed in ARPES 
from epitaxial graphene systemP^^. 

Perhaps the most interesting feature of the SL WS 2 
valence band structure is the large spin-splitting at the 
K points. This lifting of the spin degeneracy is a charac¬ 
teristic feature of the SL and due to the removal of the 
inversion sy mmetry p resent in the 2H bulk material and 
in the bilayei^^^^^®^. However, even in the SL material, 
the splitting is constrained by time-reversal symmetry. 
This forbids a splitting at T and the additional combina¬ 
tion with translational symmetry also forbids a splitting 
at M. In a simple picture, the size of the splitting at K 
strongly depends on the atomic spin-orbit splitting and 
it is thus expected to be significantly larger for W than 
for Mo. 

The spin-splitting near K is investigated in more de¬ 
tail using the data in Fig.[^ which shows a spin-splitting 
of the valence band extrema of =419(11) meV 

(Fig. [^b)). This does indeed greatly exceed the value 
observed for SL MoS 2 /Au(lll) (AFy^=145(4) meVj^, 
and is in good agreement with theoretical predictions^ 
and with values measured for analogous exfoliated 
materialJ^. Among the layered TMDCs, only SL WSe 2 
(AFy^=462 meV) and SL WTe 2 (A£'y^=480 meV) are 
expected to exhibit a larger spin-splittin ^^^ l ^^ , something 
that was recently confirmed for SL WSeJ^. 

A closer look at the contours in the constant energy 
surfaces (Fig. [^c)) reveals an anisotropy of the valence 
band dispersion around visible as a trigonal warping 
(TW). This effect is theoretically predicted and in the 
simplest approximation can be described as a third order 
correction to the parabolic energy dispersiorP^. In the 
k • p formalism, it is caused by the interaction terms be¬ 
tween the np perrn ost valence band with the lower lying 
valence band^^^^. TW reflects the underlying three-fold 
rotational symmetry of the crystal structure. It should 
be noted that this effect is connected to the general elec¬ 
tronic structure rather than to the relativistic spin-orbit 
coupling: thus, for example, it is also present for gr aphite 
and graphene where spin-orbit coupling is negligible^^l^. 

Given the constraint that the spin-splitting has to van¬ 
ish at T and M, the size of the splitting at K can 
be expected to indirectly affect the band curvature of 
the valence band maximum and thus the hole effective 
mass of the material. This is investigated by determin¬ 
ing the effective masses near K. For consistency with 
calculation^^ and to avoid the problem of a directional 
dependence due to the TW, the effective mass has been 
fitted in a region very close to the K point 0.07 A“^) 
along the f-K-M direction. This procedure leads to hole 
effective masses for VBl and VB2 (defined in Fig. i a)) as 


"^VBi = 0.40(02)me and mvB2 = 0.57(09)me. The same 
fitting procedure applied for SL MoS 2 /Au(lll]P^ yields 
t^vbi = 0.55(03)me and 771^^2 = 0.67(04)me. The effec¬ 
tive mass of the uppermost valence band is therefore in¬ 
deed reduced, consistent with the naive expectation (even 
though this model does not attempt to explain why mvB 2 
is also reduced). All these values are in good agreement 
with calculation^^. 

Finally, an interesting difference between the data 
presented here and that reported previously for SL 
moS2/Au(iii;I^ is the very clear presence of the 
Au(lll) surface state near T seen in Fig. For SL 
MoS 2 /Au( 111 ), only a very faint signature of this state 
was reported. This difference is ascribed to the slightly 
different preparation procedure. In the present work, the 
synthesis and all the analysis were performed in a single 
UHV system without ever exposing the sample to air. In 
the previous work for SL MoS 2 /Au(lll), synthesis and 
ARPES were performed in separate UHV systems, and 
the sample was transferred between them through air. 
After such a transfer, atomically clean SL M 0 S 2 can be 
recovered by a brief anneal, and the measured electronic 
structure of the layer is not affected. This, however, is 
not necessarily the case for the remaining clean Au(lll) 
terraces, where the surface state might remain quenched 
by adsorbed contaminants — at least, this may be the 
case for the low annealing temperature of 500 K used in 
Ref. [25] The observed surface state in Fig. is thus 
likely to be located on the remaining clean terraces and 
not under the SL WS 2 . 


IV. CONCLUSIONS 

We have introduced a preparation method for high- 
quality epitaxial SL WS 2 on Au(lll) and studied the 
electronic structure of this system. The observed va¬ 
lence band dispersion of SL WS 2 /Au(lll) around K 
and the spin-splitting at K of AFy^ =419(11) meV are 
found to be in good agreement with calculations for the 
free-standing layei^. The strong spin-orbit splitting con¬ 
tributes to a lowering of the hole effective mass near the 
valence band maximum, suggesting that WS 2 could be 
a more suitable material for electronic applications than 
M 0 S 2 . 
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(a) hv=45 eV (b) hv=82 eV 
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FIG. 4: (color online) (a) ARPES spectra of SL WS 2 on Au(lll) along high-symmetry directions shown in the inset. The 
well-defined contribution from the An surface state (SS) and the sp bands are indicated in the image. The projected bulk band 
gaps of Au(lllj^are outlined by dashed black lines. The energy differences between the valence band local extrema at the f 
and K (0.51 eV) points and spin-splitting of the bands at K (0.42 eV) are marked. The theoretical dispersion for free-standing 
SL (red dashed line) is superimposed on the data (after Ref. fTS]) . (b) Photoemission intensity along constant energy contours 
throughout the surface Brillouin zone at different binding energies. 



FIG. 5: (color online) (a)ARPES spectrum of SL WS 2 on 
Au(lll) along the A-K-A direction, as indicated in the in¬ 
set, (b) Energy distribution curve through the K point (red 
dashed line in (a)) fitted with two Lorentzian profiles and 
quadratic background, (c) Gonstant energy contours close to 
the K point. The Qx and qy coordinates have K as origin. 
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